UdeM-LPN-TH-93-156 
NYU-TH-93/05/02 
July 1993 



On the asymptotic 0{aas) behavior of the electroweak gauge bosons 
vacuum polarization functions for arbitrary quark masses. 

A. Djouadi^I and p. Gambino^. 

^ Laboratoire de Physique Nucleaire, Universite de Montreal, Case 6128 Sue. A, 

H3C 3J7 Montreal PQ, Canada. 

^ Department of Physics, New York University, 4 Washington Place, 
New York, NY 10003, USA. 



ABSTRACT 

We derive the QCD corrections to the electroweak gauge bosons vacuum polar- 
ization functions at high and zero-momentum transfer in the case of arbitrary 
internal quark masses. We then discuss in this general case [i) the connection 
between the 0{aas) calculations of the vector bosons self-energies using di- 
mensional regularization and the one performed via a dispersive approach and 
{ii) the QCD corrections to the p parameter for a heavy quark isodoublet. 
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1. Introduction 



The Minimal Standard Model of the electromagnetic, weak and strong interactions [with 
three generations of fermions and one Higgs doublet] has achieved a tremendous success in 
describing all experimental data within the range of energies available today. In particular, 
it has been tested to the level of its quantum corrections in the high-precision LEP and SLC 
experiments |I| . To allow for such precision tests, a considerable amount of theoretical work 
has been done in the last few years to calculate the relevant radiative corrections 0. 

Among these radiative corrections, a special attention has been devoted to the effects of 
the top-bottom weak isodoublet in the vacuum polarization functions of the electroweak 
gauge bosons. Virtual top quarks appear in these self-energies and lead to contributions 
that are quadratically proportional to the large mass of this still undiscovered particle 0]. 
These leading quadratic corrections can be parametrized by a deviation from unity of the 
p parameter |^ which measures the relative strength of the neutral to charged currents at 
zero-momentum transfer. 

In fact, due to the large size of these contributions and in view of the remarkably high 
accuracy of the experimental data, some corrections were required beyond the one-loop 
approximation. The first one involves virtual Higgs boson exchange which leads to potentially 
large Oi^a^m^ / M^r) contributions ||^, |^. Motivated by the large value of the QCD coupling 
constant, the second type of corrections are the two-loop 0{aas) contributions which have 
been calculated in using dimensional regularization and in |^ via a dispersive approach; 
the connection between these two calculations [which is needed in order to implement non- 
perturbative effects in various renormalization schemes] has been discussed recently in [^, p!0|] . 

Once the top quark is found, a comparison of the experimental data with the accurate 
theoretical predictions will allow to search for virtual effects of New Physics beyond the 
Minimal Standard Model. An immediate candidate would be a fourth generation of fermions, 
the existence of which is still allowed by present data if the associated neutrino is heavy 
enough [|l|. The latter would contribute to the vector bosons self-energies in a way similar 
to the t-b isodoublet. For instance, if the mass splitting between the two quarks is large, 
the main contribution will be cast into a correction to the p parameter. 

However, the previously mentioned two-loop results cannot be extrapolated to this case: due 
to the technical difficulties that one encounters when calculating at this level of perturbation 
theory, these results involve the simplifying approximation that the bottom quark is massless. 
[Note that this is not the case of the Higgs correction to the p parameter derived in |^ where 
the full dependence on has been retained, but the Higgs boson mass has been neglected, 
an approximation which turned out to fail badly [0 for realistic values of the quark masses] . 
While this approximation is certainly good enough for the top-bottom isoudoublet, it is very 
unlikely that this will be the case for extra quark isodoublets. 
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In this paper, we calculate using dimensional regularization the two-loop QCD corrections to 
the electroweak gauge bosons vacuum polarization functions at high-momentum and at zero- 
momentum transfer in the case of arbitrary quark masses. This allows to derive the 0{aas) 
corrections to the p parameter, to which heavy quark isodoublets will mainly contribute 
if there is a sizeable splitting between the two quark masses. The full QCD corrections, 
including the subleading terms, can be obtained via a [numerical] dispersive integral of the 
imaginary parts of the vector boson self-energies which are known in the case of arbitrary 
quark masses [O]. Using the high-momentum transfer behavior that we derive here, one 



can obtain the vacuum polarization functions in dimensional regularization; this allows to 
implement non-perturbative effects in any renormalization scheme. 

The paper is organized as follows. In the next section we summarize the one-loop results 
which will be relevant to our discussion. In section 3, we discuss the high-energy behavior of 
the vacuum polarization functions and generalize to the case of arbitrary quark masses and at 
the two-loop level, recent analyses of the connection between calculations using dimensional 
regularization and the one using dispersive integration p. In section 4, we discuss the 
behavior of the vacuum polarization functions at low-momentum transfer and derive the 
QCD corrections to the p parameter in the case of a heavy quark isodoublet. Section 5 
contains our conclusions. 



2. One— loop results 



To set the notation and for the sake of completeness, we summarize in this section all the 
one-loop results which will be relevant to our discussion. 

The contribution of a fermionic loop to the vacuum polarization tensor of a vector boson i, 
or to the mixing amplitude of two bosons i and j, denoted 11^^, is defined as 



(fxe"'-'' < OIT* 



10 > 



(2.1) 



where T* is the covariant time ordering product and q the four-momentum transfer; J^, 
are fermionic currents coupled to the vector bosons i,j and constructed with spinor fields 
whose corresponding masses are ma, mb- The vacuum polarization tensor can be decomposed 
into a transverse and a longitudinal part. 



n?(g^) + 



(2.2) 



and the two components can be directly extracted by contracting 11^^ (g^) by the Lorentz 
invariants g^^ — q^qu/'f and q^qu/q^] they can be written as [with s = g^] 

vV n^,i(s) + a'a^ n^,L(s) ] (2.3) 



n^.L( 



a 

IT 
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with f * and a* the vector and axial-vector couphngs of the gauge boson i to the fermions 
expressed in units of the proton charge e = yAna. Note that, for amphtudes involving very 
light external fermions, only the transverse part of the self-energies contribute significantly. 

The vacuum polarization function can be calculated directly by evaluating, for instance at 
the one-loop level, the diagram of Fig. la using dimensional regularization. Alternatively, 
the polarization function can be expressed as a dispersive integral of its imaginary part 
[which is related to the decay widths of the vector bosons into fermions], 

ns,A =- / ^ 2.4 

TT J{ma+mt)^ S — S' — 16 

where the cut-off scale A 3> y/s,m,a,mh regulates the ultraviolet divergence. However, in 
order to be meaningful, the previous dispersion relation needs to be subtracted. For a 
conserved current, as in the case of QED, the presence of spurious quadratic divergence is 
removed by a subtraction of the integral at s = 0, restoring in this way the Ward identity 
which requires the photon vacuum polarization function to be transverse and proportional 
to s. The situation becomes more complicated for a non-conserved current, for which some 
logarithmic divergencies are expected to be proportional to the squared fermion masses, and 
should not be removed in the subtraction procedure. 

Recently, Kniehl and Sirlin proposed a subtraction prescription for the dispersion relation 
eq. (2.4) which manifestly preserves the generic Ward identities constraining the vacuum 
polarization functions Let us recall the basic elements of this procedure and write the 
Ward identity 

q,A{s) = I " < 0|T* [a^J^(x) 4(0)] |0 >= g'^n^.(s) (2.5) 

with a notation as in eq. (2.1) but where we dropped the indices i,j for simplicity. Contract- 
ing the previous equation by g^, one obtains A(s) = IIl{s)] then using the linear combination 
s$(s) = n^(s) — nj'(s), one can rewrite the transverse part in terms of A(s) and $(s), 

Ut{s) = -s^s) + A{s) , A{s) = Ul{s) (2.6) 

Now $(s) and A(s) are only logarityhmically divergent for infinite energies, unlike nT('S) 
which is quadratically divergent as discussed above. For $(s) this is obvious for dimensional 
reasons [two powers of the external momentum have been extracted], and for A(s) this is 
because the soft breaking of the current by mass terms implies that its derivative involves 
at most operators of canonical dimension three, and one power of the momentum is extracted. 

Following the authors of Ref. we therefore write the once subtracted dispersion relation for 
the transverse part of the vacuum polarization function expressing nr(s) through eq. (p.6|) , 
and then writing unsubtracted dispersion relations for nx(s) and $(s). 



Ut{s) = ds' ^"^^^Tl^'J +Xm^{s') +— / '^mL{s')-s^s')] (2.7) 

^' nJim.+m^r s'-s-ie ^ ' 27rz/|s'i=A2 s' ^ ^ ' v ;j v ; 



ImI\.T{s' 



1 r ds' 

|s'|=A2 
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where the quadratic divergence is now cancelled by the contribution of Xm$(s'), and the 
logarithmic divergence is regulated by the cut-off A. The integral over the large circle, where 
the s in the denominator has been neglected compared to s', comes from the fact that 
and do not vanish asymptotically but behave like constants modulo logarithms. Note 
that for conserved currents d^J^{x) = 0, IIl{s) vanishes reducing then the Ward identity to 
the QED expression IIt{s) = — 

One can now express the result of the dispersion integral eq. (2.7) for arbitrary quark masses 
and for vector and axial-vector currents, as 

U^'^is) = sX, + {ma T m,fY, + F^'\s) + ^ ^ [t^V{s') - s^^'^^s')] (2.8) 

where the function Fi''^{s) behaves like a constant as s — > oo. Xi and Yi are divergent 
constants involving logarithms of the cut-off scale A; their explicit form is given by 

Xi = log ^, ri = -^log + 3 (2.9) 

Analogously, the real part of H^'^^i^s) calculated in n = 4 — 2e dimensions can be written as 

n^'^(s) = sX, + {rua T mt^Y^ + F^'^{s) (2.10) 

where now, the divergent constants Xi and Yi involve poles ine = 2 — n/2 and logarithms 
of the 't Hooft mass scale fi; their expressions at the one-loop level are given by 



~ 1 maTTib ~ 3 3 rriamb , . 

X, = --log^, >^i = -^+2log^ (2.11) 



The finite function F^''^{s) in eq. (2.8) is, of course, the same as in eq. (2.10). Since the result 



of eq. ( |2.7| ) should be consistent with the one obtained using dimensional regularization, the 
integral over the large circle must be absorbed by a redefinition of the divergent constants 
Xi and Yi. In particular, if IIl{s) and ^(s) are evaluated in dimensional regularization, the 
effect of the second integral in eq. (|2.7|) should be simply to replace Xi, Yi by Xi, Yi. This 
is indeed true and one has 

1 r ^^'^(s') 
^{ ds'^^ = X,-X, 

2m 7|s'|=A2 s' 



ds' ^ ; ^ = (m,Tm,)^(ri-Fi) (2.12) 
27n J|s'|=A2 s' 

This point is crucial to ensure that the cancellation of divergences occurring in physical 
observables like the p parameter yields a perfect independence of these observables from the 
calculational approach used to obtain the vacuum polarization functions. 
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Before we generalize this analysis for arbitrary quark masses at the two-loop level, let us 
briefly discuss the contribution of a heavy quark isodoublet to the p parameter. 

The fermionic contribution to the p parameter can be expressed in terms of the difference 
between the transverse components of the Z and W bosons self-energies at zero-momentum 
transfer, = 



i-1 

P 



(0) 



(0) 



(2.13) 



Besides the light fermions contribution to the photon vacuum polarization function which 
can be mapped into the running of the fine structure constant a, Ap represents the leading 
part of the radiative correction to electroweak observables such as the W boson mass and 
the effective weak mixing angle sin^ 6w as measured for instance on the Z resonance . 



The contribution of a quark pair with different masses, nia 7^ nib 7^ 0, to 11^^(0) in the 
one-loop approximation. Fig. la, is 



nK(o) 



[m. 



+ 1- Pa- Pb] + 



4 I 4 

rri^. + mi 



m^a-^b 



log 



mi 



± Smamf, 



1 

- + 1 

e 



-{Pa + Pb) - - 



1 mi + m? 



2 — ml 



1 



mt 



(2.14) 



with pafi = log m^f,//i^ where p is the 't Hooft mass scale introduced to make the coupling 
constants dimensionless in n = 4 — 2e dimensions [note that we have introduced an extra term 
{e'~'/4:7ry, where 7 is the Euler constant, to prevent uninteresting combinations of log47r, 7, • • • 
in the previous result]. One can easily verify that for a vector current with equal fermion 
masses = m?,, as it is the case for the photon or the 7-Z mixing amplitude, the vacuum 
polarization function vanishes at = 0. 

In the case of a heavy quark isodoublet [with Standard Model couplings] Ap can alternatively 
be written in terms of the vector and axial-vector components of the vacuum polarization 
function eq. (2.3); using the fact that 11^(0, ma, ma) vanishes and that in this limit the 
transverse and longitudinal components are equal, one has 

V2Gf 



Ad) 



P 



-fi{s = 0,ma,mb) 



(2.15) 



f{s,ma,mb) 



1 



(s, ma, ma) + n^(s, mh, m,,) - n^(s, ma, mi,) - n^(s, ma, m^) 



Using the previous expressions of 11^'^ (0) one obtains the well-known one-loop expression 
that is quadratic in the fermion masses [HI 



/i(0, ma, mfe) 



ml + ml + 



2mlml ml 
log 



mi 



mt 



mi 



(2.16) 
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3. High— momentum transfer 

As previously discussed, one can calculate the QCD corrections to the vacuum polarization 
functions of the electroweak gauge bosons using a dispersive approach: the polarization 
functions can be expressed as dispersive integrals of their imaginary parts, the expressions of 
which are available in the literature at this order of perturbation theory for arbitrary quark 
masses [0|. However, to connect this result with the full expression that can be derived 
using dimensional regularization, one needs the high-momentum transfer behavior of the 
polarization function in the latter scheme. 



We have therefore calculated the expressions of the vacuum polarization functions of the 
weak gauge bosons at order 0{aas), Fig. lb plus the corresponding counterterms, in the 
limit of high-momentum transfer [and also at zero-momentum transfer as will be discussed 
later] for arbitrary quark masses 7^ mf, 7^ in the on-shell mass scheme. This scheme, 
where the quark masses are defined as the poles of their propagators, is usually used to 
calculate electroweak radiative corrections |12[. Here we will simply give the final result; 
details on the calculation will be found in Ref. |jl3| . 

Using the same notation for the vacuum polarization functions as in the one-loop case, and 
defining the latter similarly to eq. (2.3) 



a as 

IT TT 



vVU^As) + a'am^r{s) 



(3.1) 



one obtains for n^'^(s) in the limit |s| ^ 00 



H^'^^dsl 00) 



nl'^d^l - 00) 



+ 
+ 
+ 



+ 



11 1 55 

^ - 2 (P« + Pb) + ^{^oga + log/3) + - - 4C(3) 



{nia =F mbf 



3 1 

2^^ Ye 
9 



11 

y 



3pa - 3p; 



11 
T 



{pa + pb) 



11 



{pa + Pfe) - 7 (log a + log /3) - - log a log /3 + 6C(3) 



[m„ 



2\ 1 ^a 

ml) log 

mi 



3 3 3 3 

2(P« + Pfc) -4(log« + log^) + 4 



+ 3(m^ + m^)(loga + log/3) 



[ma =F mi,f 
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11 , ,3 

-^{Pa + Pb)-- 



+ -Apa + Pb)-- (log a + log /?) - - log a log /3 + 6C(3) 



m 



mt 



2 r 



-^^ + -^{Pa + Pb) 



- (log a + log /3) - 5 



(3.2) 



(3.3) 



with a = —m^^/s, (3 = —ml/s and ({3) = 1.202. At this stage a few remarks are mandatory. 
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i) In the previous expressions the momentum transfer is defined in the space-hke region, 
s < 0. The continuation to the physical region can be straightforwardly obtained by 
adding a small imaginary part —ie to the quark masses squared. This reduces to make 
the substitutions log(m^ {,/ — s) ^ log — s| + in. 

a) As expected, only the transverse part of the vacuum polarization function is quadrat- 
ically divergent for |g| oo. This divergent term, the expression of which is in 



agreement with the one obtained in Ref. [^, is the same for axial and axial-vector 



currents as expected from chiral symmetry. Moreover, as required by the Kinoshita- 
Lee-Nauenberg theorem [|14| , this term does not introduce any mass singularity as nia^b 
tend to zero: — (pa + Pb) and log a + log/5 combine to give 21og(/i^/ — s). 

in) Since the vector part of the longitudinal component should vanish for = rrih, it must 
be proportional to {ma — rrih) or (m^ — ml); and since the axial-vector component can 
be obtained by changing the sign of one of the two masses, it must be proportional 
{ma + mb) or (m^ — ml) [log m^/m^ alone would have introduced mass singularities]. 
This behavior is explicitely exhibited by the previous expression of 11^'^ {s). 



Let us now write the real parts of n^''^(s), analogously to eq. (2.10), 

n^'^(.) = sX2 + {ma T m,)^% + {ml - ml) log Z, + F^'\s) (3.4) 

mi 

where the two-loop divergent constants X2, I2 and Z2 involve only the poles in e and the 
logarithms of the scale /x of the expressions of eq. (3.3); they are given by 

~ 1 mamb 

X2 = —-log — 

2e fi^ 

~ 3 1 /ll mamb\ 11, mamb , 2 "^^"^6 
Y2 = + - — - 3 log — - — log — + 3 log 



2e2 ' e V 4 P^ J 2 ^ jj? 

— + 3 log 

2e /i^ 



Z2 = -— + 3 log ^ (3.5) 



Using the results of the imaginary parts of the vector bosons vacuum polarization functions 
at 0{aas) given in Ref. |TT|] for arbitrary quark masses, one can express the result of the 
dispersive integral eq. (2.7) in this case as 



n^'^(s) = sX2 + {ma T mb)'Y2 + {ml- ml) \og^ Z2 + F^'^{s) 

mi 



2lTi J\s'\=A^ s' 



(3.6) 
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where the finite function^ F2'^{s) should be the same as the one in eq. (3.4). This has 
been verified up to 0{aas) @] in the special cases for which both approaches were used, i.e. 
'^a = and rrih = 0. 

Here again, since the result of eq. (3.6) should be consistent with the one obtained in dimen- 
sional regularization, the integral over the large circle must be absorbable in a redefinition 
of the divergent constants ^2,1^2 and Z2. Indeed, when we use dimensionally regularized 
expressions for <l>(s) and IIl{s), the integral involving $(s) in eq. (3.4) will replace X2 by X2, 
and the one involving 11^(5) will replace Y2 and Z2 by respectively Y2 and Z2. One would 
have 

ds'^-^ = X2-X2 

i ds'^^ = {m,Tni,)\Y2-Y2) 
2m 7|s'|=A2 s' 

2 

+ {ml-ml)\og^{Z2-Z2) (3.7) 

Using the 0{aas) results of the imaginary parts of the vacuum polarization functions given 
in and after a straightforward computation, we obtain for the two-loop constants X2, 12 
and Z2 

X2 = log + 4C(3)-- 

mamb 12 

3, 2 A2 9, A2 3, 3 , vr^ 

Y2 = -log' -log iog2^+ +6C(3)-- 

2 mamb 2 mamb 2 mb o 2 

Z2 = -^log^^ + 5 (3.8) 

2 mamb 

The divergent part [as A tends to infinity] of these constants for arbitrary masses, as well as 
their expressions in the special cases = mb and = have been derived in Ref. §]. 

As a consequence of the use of eq. (3.7) in the dispersion relation eq. (2.7), one can now 
implement in the general case of arbitrary quark masses, non-perturbative effects |P, |10[ in 
the absorptive part of the amplitude [that can be calculated by integrating, e.g. numerically, 
the expressions of the imaginary parts of the polarization functions derived in Ref. [pJj 



defined not only in the on-shell scheme, which is natural in the dispersive approach, but 
also in the MS scheme for which a dimensional regularization calculation is required. 



^Note that the definition of F^'^^is) used here differs from the one used in Ref. since we normahzc the 
divergent constants to the dimensional regularization result eq. (3.2) keeping only the ultraviolet divergencies 
and the related logarithms. 
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4. Zero— momentum transfer 

We come now to the discussion of the zero-momentum transfer behavior of the vacuum 
polarization functions. In this hmit, the evaluation of the two-loop diagrams Fig. lb in the 
on-shell mass scheme, leads to the following expression for n^^^(O) 



n^^(o) 



3,9 2\ 



^{ml + ml) - 3mlpa - Smlpb 



+ 3^aPa + 3m;,P(, + — (m„ + 7716) + - (m^-mft) 



+ 
± 



G 



2 2 
mtm 



mi 



log3| + m>g " ' 



mt 



mamb 



mi 
3 1 

■- + - 

e 



log'^ 
(m2 - mlY ml 



3pa + 3pfe - 



llA 3, ,2 



+ 3 



+ 3- 



2 2 
mtmt 



m„ 



mi 



m^)2 



mt 



mt 



mi 



4(Pa + Pb) - 



31 



(4.1) 



where, in terms of the Spence function defined by Li2(a:) = — lo V ^log(l — xy)dy, the 
function G is given by 



G{x) = 2U2{x) + 21ogxlog(l -x) + 



X 



1 — X 



log^ X 



(4.2) 



As one might have expected, there is no singularity in the self-energies in this limit. In 

addition, besides the manifest symmetry in the exchange ^ mt, the previous expres- 
sions exhibit the facts that n^'^(O) can be obtained from ll^']^{0) by simply making the 
substitution ma{mb) —>■ —ma{—mb) as expected from 75 reflection symmetry and that the 
longitudinal and transverse components are equal. These features provide good checks of 
the calculation. 

Using the previous expression, one readily obtains the QCD corrections to the contribution 
of a heavy quark isodoublet to the p parameter in the general case m^ ^ mb ^ 0. Defining 
the contribution to the p parameter analogously to eq. (2.15) 



A(^V= ^^/2(0, m„,m,) 



87r2 



TT 



(4.3) 



the function /2 will be given by 



f2{0,ma,mb) 



-3 < + 2 ^2 



mlml 



X 



a ^b 
2 



log 



m: 



mt 



1 + 



ml + ml 



mt 



2U2 



'ml\ ml^ / mV 

+ 21og-^log 1- 



mt 



m„ 



log 



mt 



mt 



mt 



ml — ml 



log 



mg 
ml 

2 



+ (m„ - m^) 



mt 



y 



(4.4) 
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/2(0,ma,mb) is free of ultraviolet divergences as it should be since Ap is an observable 
physical quantity; furthermore it does not depend on the 't Hooft mass scale fi. Note that 
the symmetry in the interchange of rria and mi, is now hidden in the term in the last line 
of the previous equation, but this term is simply G(m^/m^) — G(m^/m^) and we have used 
the fact that G{l/x) + G{x) = 27rV3. 

In the limit of large mass splitting between the two quarks, rua S> mb, the QCD corrections 
to the p parameter reduce to the known result for = 



The contribution of a new quark isodoublet to the p parameter is exemplified in Fig. 2a 
both in the one-loop approximation [full line] and including the QCD radiative corrections 
[dashed line]. We have fixed the mass of one quark to 250 GeV and varied the mass of 
the other quark from 100 to 400 GeV; for the strong coupling constant we have used the 
numerical value as — 0.12. [The scale at which as should be evaluated is the mass of the 
heavy quark |jl6|; this choice can also be justified by arguments based on effective field theory 



As one can see the contribution of the heavy quark isodoublet to Ap is very small for 
approximately mass degenerate quarks, ma = mi,] this illustrates the well known fact that 
large contributions to Ap are induced only if the splitting between the two quark masses is 
large. We see that the QCD corrections are always negative and therefore reduce the size of 
Ap. In fact, in this mass range the QCD corrections are practically constant and amount to 
a decrease of Ap by approximately 15% as can be seen from Fig. 2b. 

To close this section, let us briefly discuss the application of the dispersive method to physical 
observables, taking the example of the fermionic contribution to the p parameter. Using the 
fact that at zero-momentum transfer, the transverse and longitudinal components of the W 
and Z self-energies are equal, and recalling that the vector component vanishes in this limit 
for equal mass quarks, we have written Ap in terms of the function f{s,ma,mi,) defined 
in eq. (2.5). Since Ap is an observable physical quantity [and therefore free of ultraviolet 
divergences], and since the divergences of 11^ are independent of s, f{s,ma,mi,) can be 
evaluated by means of an unsubtracted dispersion relation 



We have verified at the two-loop level, that the ultraviolet finite function /(s, m^, m^) tends 
to zero as s becomes infinite and therefore satisfies the condition 




(4.5) 




(4.6) 




(4.7) 



which generalizes the result obtained in Ref. |^ in the case of one massless quark, mb = 0. 
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5. Summary 



In this paper, wc have derived the asymptotic behavior of the heavy quark contribution to 
the electroweak gauge bosons vacuum polarization functions at 0{aas) for arbitrary quark 
masses, nia 7^ 7^ 0. 

Using the high-momentum transfer behaviour, we have extended to the case of arbitrary 
fermion masses at the two-loop level, recent analyses of the connection between calculations 
using dimensional regularization and the ones performed by means of dispersive integration. 
One can therefore implement, in the general case, non-perturbative effects in the absorptive 
part of the vacuum polarization function [which can be calculated using the available ex- 
pressions of the imaginary parts] not only in the on-shell scheme but also in the MS scheme 
which inherently requires a dimensional regularization calculation. 

With the results that we obtain for the zero-momentum transfer behavior of the electroweak 
gauge bosons self-energies, we derived the QCD corrections to the leading contribution of 
new heavy quarks to electroweak observables: the deviation of the p parameter from unity. 
This correction is practically constant and amounts to a decrease of the one-loop value of 
Ap by approximately 13-15% for as ~ 0.12. 
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Figure Captions 

Fig.l Contribution of a quark pair to the vacuum polarization function of a gauge boson in 
the one-loop approximation (la) and at the two-loop level (lb). 

Fig. 2 Contribution of a heavy quark isodoublet to the p parameter. (2a) At 0{a) [full line] 
and total contribution at 0{aas) [dashed line], (2b) total contribution at 0{aas) 
normahzed to the one at 0{a). We have taken as ^ 0.12. 



14 



